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CAMENEN,  B.  and  LARSON,  M.,  2008.  A  general  formula  for  noncohesive  suspended  sediment  transport.  Journal  of 
Coastal  Research,  24(3),  615-627.  West  Palm  Beach  (Florida),  ISSN  0749-0208. 

A  simple  and  robust  suspended  load  transport  formula  for  noncohesive  sediment  is  presented  for  application  to  river, 
estuarine,  and  coastal  environments  with  the  use  of  depth-averaged  models.  The  formula  is  based  on  an  exponential 
profile  for  the  concentration  and  assumes  a  constant  velocity  over  depth  to  simplify  the  calculations.  These  assump¬ 
tions  were  validated  with  a  large  data  set,  including  data  with  a  steady  current,  wave  and  current  interaction,  and 
breaking  waves.  The  formula  has  two  parameters:  the  mean  sediment  diffusivity  over  depth  and  the  bottom  reference 
concentration.  The  sediment  diffusivity  is  estimated  assuming  a  linear  combination  of  mixing  because  of  breaking 
waves  and  the  energy  dissipation  in  the  bottom  boundary  layer  from  the  mean  current,  waves,  or  both.  The  bottom 
reference  concentration  is  a  function  of  the  Shields  parameter.  Overall,  the  formula  developed  in  this  study  yields  the 
best  agreement  with  the  compiled  data  set  compared  with  a  number  of  existing  formulas  for  estimating  the  suspended 
load. 


ADDITIONAL  INDEX  WORDS:  Suspended  load,  concentration,  profile,  reference  concentration,  diffusion  parameter, 
current,  waves,  breaking  waves,  noncohesive  sediment. 


INTRODUCTION 

Accurate  prediction  of  noncohesive  sediment  transport 
rates  is  an  important  element  in  morphological  studies  for 
river,  estuarine,  and  coastal  environments.  Depth-averaged 
(2DH)  models  are  widely  employed  nowadays  and  generally 
allow  for  the  use  of  one  sediment  transport  formula  only  in 
the  entire  domain.  Such  a  formula  needs  to  be  robust  and 
reliable  in  a  wide  range  of  conditions.  In  estuarine  and  coast¬ 
al  environments,  the  process  of  sediment  transport  becomes 
increasingly  complex  because  of  the  presence  of  oscillatory 
flows  and  the  interaction  between  steady  and  oscillatory 
flows.  For  example,  for  longshore  sediment  transport,  the  in¬ 
fluence  of  short  waves  is  expressed  as  sediment  stirring, 
which  increases  bed  shear  stress  and  the  vertical  mixing  co¬ 
efficient  (diffusion)  for  sediment  in  suspension  (Bijker,  1968; 
Van  Rijn,  1993;  Watanabe,  1982).  The  development  of  prac¬ 
tical  sediment  transport  models  still  has  a  strong  empirical 
character  and  relies  heavily  on  physical  insights  in  combi¬ 
nation  with  quantitative  data  obtained  in  laboratory  and  field 
studies. 

The  earliest  formulas  were  mainly  based  on  the  concept 
that  the  sediment  transport  rate  for  steady  uniform  flows  can 
be  related  to  bottom  shear  stress  (Einstein,  1950;  Enge- 
lund  and  Hansen,  1972;  Meyer-Peter  and  Muller, 
1948)  and  assumed  that  bed  load  transport  prevailed.  How- 
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ever,  when  the  bottom  shear  stress  is  large  enough,  the  sed¬ 
iment  particles  can  be  lifted,  put  in  suspension,  and  trans¬ 
ported  in  large  quantities  by  the  current.  Thus,  suspended 
load  is  often  dominant  for  fine  sediment  (median  grain  size 
d50  <  0.5  mm)  under  medium  shear  stress,  under  the  pres¬ 
ence  of  bed  forms,  or  with  wave  stirring.  The  depth-averaged 
volumetric  suspended  load  transport  qsg  is  herein  defined  (see 
Nielsen,  1992,  p.  201)  as  the  integrated  product  of  the  ve¬ 
locity  u  and  the  concentration  c  from  the  edge  of  the  bed  load 
layer  (z  =  zR)  to  the  water  surface  (z  =  h),  averaged  in  time, 
yielding, 


fh 

qss  =  c(z,  t)u(z,  t)  dz  (1) 

JzR 

where  h  is  the  water  depth,  z  is  a  vertical  coordinate,  and  x 
is  the  time-averaged  value  of  the  variable  x. 

Assuming  that  variables  u  and  c  can  be  decomposed  into 
two  components,  i.e.,  a  time-averaged  component  u  and  c  and 
an  oscillating  component  u  and  c,  u(z,  t)  =  u(z)  +  u(z,  t) 
and  c(z,  t)  =  c(z)  +  c(z,  t),  then  Equation  (1)  becomes: 


fh  fh 

qss  —  c(z)m(z)  dz  +  c(z,  t)u(z,  t)  dz  (2) 

JzR  JzR 

where  the  first  term  corresponds  to  the  current-related  sus¬ 
pended  load  and  the  second  term  to  the  wave-related  sus- 
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Figure  1.  Typical  concentration  profile  (full  line)  from  the  fixed  bed  to 
the  water  surface  on  the  basis  of  experimental  results  by  Dohmen-  Jans¬ 
sen  and  Hanes  (2002),  together  with  the  exponential  (dashed  line)  and 
power  law  (dash-dotted  line)  profiles  fitted  to  the  suspended  concentra¬ 
tion. 


pended  load.  Steady  conditions  are  generally  assumed  to  sim¬ 
plify  the  problem,  so  time-averaged  values  u(z)  and  c(z)  are 
used.  Therefore,  an  accurate  estimate  of  the  total  suspended 
load  requires  correct  predictions  of  the  mean  current  velocity 
and  the  concentration  profile.  Even  if  the  unsteady  part  of 
the  suspended  sediment  transport  is  significant  for  some  cas¬ 
es  (phase  lag  effects  over  ripples;  see  Van  der  Werf  and 
Ribberink,  2004),  in  this  study,  we  will  focus  on  the  steady 
part  because  most  of  the  experimental  data  measured  this 
quantity  only. 

The  overall  objective  of  this  study  was  to  develop  a  reliable 
and  general  formula  for  predicting  suspended  load  transport 
under  a  wide  range  of  river,  estuarine,  and  coastal  conditions. 
Such  a  formula  is  proposed  for  2DH  models  in  which  only  an 
average  value  of  the  velocity  over  the  water  depth  is  provid¬ 
ed.  This  paper  is  organized  as  follows:  The  proposed  formula 
is  first  developed,  then  the  experimental  data  used  for  this 
study  are  described,  together  with  a  validation  of  the  hypoth¬ 
esis  underlying  the  proposed  formula.  The  study  and  calibra¬ 
tion  of  the  two  main  parameters,  which  are  the  reference  con¬ 
centration  cR  and  mean  diffusivity  e,  are  presented.  Finally, 
conclusions  concerning  the  derived  net  suspended  load  trans¬ 
port  formula  and  its  predictive  capability  are  discussed. 

DEVELOPMENT  OF  A  SUSPENDED  LOAD 
FORMULA 

On  the  basis  of  experimental  results  by  Dohmen-Janssen 
and  Hanes  (2002),  Figure  1  presents  a  typical  concentration 
profile  in  which  sheet  flow  and  suspended  load  coexist.  They 
observed  a  drop  in  concentration  at  the  top  of  the  sheet  flow 
layer  and  defined  the  edge  of  the  sheet  flow  (bed  load)  layer 
in  which  the  concentration  c  =  cb  =  0.08  (the  maximum  vol¬ 
ume  concentration  is  defined  as  c0  =  0.65).  It  appears  that 
close  to  the  bed  load,  both  power  law  and  an  exponential  pro¬ 


file  underestimate  the  concentration.  On  the  other  hand, 
these  two  theoretical  profiles  yield  very  similar  results  for  the 
suspended  concentration.  Indeed,  compared  with  the  experi¬ 
mental  data  compiled  (see  Experimental  Data  section),  both 
exponential  or  power  law  profiles  could  be  fitted  to  the  ex¬ 
perimental  data  with  a  small  relative  error. 

Because  the  use  of  a  power  law  profile  for  the  sediment 
concentration  requires  a  reference  level  (z  =  zR),  which  in¬ 
duces  an  additional  parameter  and  thus  even  more  uncer¬ 
tainty  (because  zR  is  often  arbitrarily  chosen),  an  exponential 
law  profile  was  preferred  (zR  =  0  could  then  be  assumed), 
assuming  a  constant  value  for  the  sediment  diffusivity  over 
depth  e.  By  solving  the  mass  conservation  equation  for  the 
steady  equilibrium  of  a  particle  under  gravity  and  hydrody¬ 
namic  forcing,  the  following  profile  for  the  sediment  concen¬ 
tration  is  obtained, 


c(z)  =  cRexp 


(3) 


where  cR  is  the  bottom  reference  concentration,  Wa  is  the  set¬ 
tling  velocity,  and  the  parameter  WJt  determines  the  sus¬ 
pension  conditions.  In  determining  qss,  following  the  simpli¬ 
fied  approach  by  Madsen,  Tajima,  and  Ebersole  (2003), 
the  vertical  variation  in  u  can  be  neglected  (also  discussed 
further  in  the  Validation  of  the  Hypothesis  section).  The  sus¬ 
pended  sediment  load  is  thus  found  to  be 

=  UcF(cR,  e)  (4) 

where  Uc  =  u  is  the  depth-averaged  velocity  and  the  function 
F  determines  the  quantity  of  sediment  available.  In  solving 
the  integral,  the  ratio  WBh/e  is  often  assumed  to  be  large, 
implying  that  the  exponential  term  exp (—WJiIe)  ~  0  or  F  ~ 

CR6/Wb. 

In  the  case  of  wave  and  current  interaction,  Equation  (4) 
can  be  modified  to  take  into  account  possible  sediment  trans¬ 
port  in  the  direction  of  the  waves  (see  Figure  2), 


1  -  exp 


=  (tfc 


cw, onshore 


-  U, 


cw, offshore- 


)F(cr,  e) 


Qssn  =  Uc  sin  9  F(cr,  e) 


(5) 


where  the  subscript  “w”  indicates  the  direction  of  the  wave 
and  the  subscript  “n”  indicates  the  perpendicular  direction; 
Ucwj  is  the  root  mean  square  value  of  the  velocity  [u(t)  =  ujjt) 
+  Uc  cos  cp]  over  the  half  period  TwJ,  for  which  the  subscript 
j  should  be  replaced  by  either  onshore  (for  u(t)  a  0)  or  offshore 
(for  u(t)  <  0);  and  ip  is  the  angle  between  the  wave  and  cur¬ 
rent  direction  (see  Figure  2).  For  a  steady  current  alone, 
Equation  (5)  reduces  to  Equation  (4),  and  for  sinusoidal 
waves  alone,  Equation  (5)  yields  zero  transport. 

Thus,  the  two  main  parameters  to  estimate  for  calculation 
of  the  suspended  load  are  the  bottom  reference  concentration 
and  mean  sediment  diffusivity  over  the  water  depth.  Here, 
various  data  sets  were  used  for  model  development  with 
steady  and  oscillatory  flows,  including  wave  breaking.  The 
aim  of  this  paper  was  to  develop  and  validate  a  formula  that 


Journal  of  Coastal  Research,  Vol.  24,  No.  3,  2008 


Suspended  Sediment  Transport  Formula 


617 


Figure  2.  Definition  sketch  for  the  depth-averaged  velocities. 


gives  accurate  results  when  both  waves  and  current  are  pres¬ 
ent  in  the  nearshore  zone.  Thus,  the  effects  of  a  steady  cur¬ 
rent,  waves,  and  breaking  waves  on  sediment  diffusivity  and 
the  reference  concentration  were  investigated  carefully  to  in¬ 
clude  all  the  main  hydrodynamic  parameters  and  provide  ro¬ 
bust  predictive  formulas. 

EXPERIMENT  DATA 
Selection  of  the  Data 

To  investigate  mean  diffusivity  over  depth,  bed  reference 
concentration,  and  resulting  suspended  sediment  transport 
in  steady  conditions,  as  well  as  for  waves  and  current  com¬ 
bined,  a  wide  range  of  existing  data  sets  were  compiled  and 
analyzed.  Depending  on  the  experiments,  velocities  were 
measured  with  impeller  flow  meters,  pitot  tubes,  acoustic 
Doppler  probes  (ADP),  laser  Doppler  velocity  meters,  or  elec¬ 
tromagnetic  current  meters.  Time-averaged  concentrations 
were  measured  with  suction  (pump)  samplers  and,  more  re¬ 
cently,  with  optical  probes  (optical  backscatter  sensors),  con¬ 
ductivity  probes,  and  ADP.  It  is  obvious  that  the  different 
instruments  and  the  precision  of  the  instruments  affect  the 
quality  and  uncertainties  of  the  results.  It  is  difficult,  how¬ 
ever,  and  not  the  purpose  here,  to  discuss  in  detail  the  un¬ 
certainties  induced  by  the  measurements  in  the  different  ex¬ 
periments.  As  pointed  out  previously,  the  limit  between  the 


bed  load  and  the  suspended  load  is  difficult  to  establish.  For 
most  of  the  cases,  concentrations  were  not  measured  close 
enough  to  the  bed,  or  the  instruments  were  not  able  to  mea¬ 
sure  very  high  concentrations.  The  suspended  load  was  esti¬ 
mated  by  the  reference  level  proposed  by  Van  Rijn  ( 1993), 
i.e.,  zR  =  max(iTr/2,  ksg),  where  Hr  is  the  ripple  height  and  kBg 
~  2 dm  is  the  grain-related  roughness  height.  An  approximate 
power  law  regression  was  used  for  the  experimental  data 
when  information  was  not  available  close  to  the  bed. 

For  a  steady  current,  Table  1  summarizes  the  data  sets  and 
lists  the  type  of  flow  motion  and  sediment  properties.  Simi¬ 
larly,  for  waves  and  current  combined,  Table  2  summarizes 
the  data  sets,  presenting  the  type  of  experiment,  wave  and 
current  conditions,  and  sediment  properties. 

For  all  the  experiments  presented  in  Table  2,  sand  with  a 
relative  density  s  =  pa/p  =  2.65  was  used  (ps  and  p  are  the 
sediment  and  water  densities,  respectively).  Most  of  these 
data  sets  were  obtained  from  the  data  compilation  provided 
by  the  SEDMOC  European  Union  research  project  (Van  Rijn 
et  al.,  2001).  For  wave-current  interaction,  only  the  cases  in 
which  the  mean  current  (preferably  the  entire  velocity  pro¬ 
file)  was  estimated  could  be  used  to  calculate  the  total  sus¬ 
pended  load.  The  number  of  measurement  points  over  the 
water  depth  are  also  shown  in  Table  2,  for  both  the  sediment 
concentration  and  the  time-averaged  velocity,  to  indicate  the 
spatial  resolution  of  the  data.  In  the  cross-shore  direction,  the 


Table  1.  Data  summary  for  experiments  on  suspended  sediment  transport  under  steady  flow. 


Authors 

Location 

Flow  Type 

No.  of 
Profiles 

d50  (mm) 

b  (m) 

Fr 

u*c  (m/s) 

Anderson  (1942) 

Enoree  River,  USA  (1940-41) 

River  data 

23 

0.7 

15 

0.15-0.25 

0.02-0.07 

Barton  and  Lin  (1955) 

Fort  Collins,  CO,  USA 

Tilting  flume 

26 

0.18 

1.2 

0.2-0. 9 

0.02-0.08 

Laursen  (1958) 

Iowa,  USA  (1961-63) 

Tilting  flume 

12 

0.4,  1 

0.9 

0.25-0.60 

0.02-0.09 

Scott  &  Stephens  (1966) 

Mississippi  River,  USA  (1961-63) 

River 

23 

0.4 

500 

0.11-0.16 

0.05-0.13 

Culberton,  Scott,  and  Bennet  (1972) 

Rio  Grande  River,  USA  (1965-66) 

River 

22 

0.18-0.33 

20 

0.2-0. 6 

0.05-0.15 

Voogt,  Van  Rijn,  and  Van  den  Berg 
(1991) 

Krammer  beach,  The  Nederlands 
(April  1987) 

Tidal  channel 

60 

0.22-0.35 

300 

0. 1-0.5 

0.03-0.15 

Peet  (1999;  see  Van  Rijn  et  al.,  2001) 

Wallingford,  UK 

Duct  experiments 

24 

0.08-0.20 

0.6 

0.2-0. 4 

0.01-0.14 

dB0  =  median  grain  size;  b  =  width  of  the  river/flume;  Fr  =  Froude  number;  u,c  =  current-related  shear  velocity. 
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Table  2.  Data  summary  for  experiments  on  suspended  sediment  transport  under  oscillatory  flow. 
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t  Values  in  the  parentheses  are  no.  of  nonbreaking  cases  and  no.  of  cases  in  which  the  suspended  load  can  be  estimated. 

$  Number  of  concentration  measurements  of  the  water  depth/number  of  time-averaged  velocity  measurements  over  the  water  depth. 
§  For  random  waves,  the  orbital  wave  velocity  Uw  is  computed  from  the  root  mean  square  wave  height  and  wave  period  Tw  =  Tp. 

**  Flat  bed  (bed  form  height  Hr  and  length  Lr  are  assumed  equal  to  zero. 
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Table  3.  Prediction  (Pred.)  of  suspended  load  transport  from  Equation  (4) 
and  the  fitted  valeu  (exponential  profile)  to  the  observed  data  for  cR  and  e 

[i(q,J  =lOg  I  %,.predl  QU.mras  V- 


Conditions 

No. 

Pred.  (%) 

X2 

X5 

Mean 

SD 

Steady  current 

187 

99 

100 

0.03 

0.09 

Waves  and  current 

322 

66 

91 

-0.16 

0.38 

Breaking  waves 

151 

72 

93 

-0.18 

0.32 

velocity  was  averaged  from  the  bottom  to  the  trough  of  the 
waves.  Thus,  it  corresponds  to  the  mean  undertow. 

An  important  question  in  nearshore  sediment  transport  is 
to  know  whether  bed  load  or  suspended  load  prevails.  For  the 
lower  regime  and  upper  regime  (sheet  flow  with  nonbreaking 
waves),  bed  load  has  been  observed  to  prevail  over  suspended 
load  (Dohmen-Janssen  and  Hanes,  2002).  For  medium  re¬ 
gimes,  when  the  bottom  is  covered  by  bed  forms  and  when 
waves  are  breaking,  it  is  often  assumed  that  suspended  load 
prevails  (Nielsen,  1992,  pp.  201-206).  With  the  Camenen 
and  Larson  (2005)  formula  for  bed  load  (not  measured  for 
all  the  collected  data),  this  assumption  was  confirmed.  For 
80%  of  the  collected  data,  the  suspended  load  largely  pre¬ 
vailed  over  bed  load  (qBJqsb  >  10),  especially  for  very  fine 
sediments  and  when  waves  were  present.  Bed  load  appeared 
to  prevail  mainly  with  coarse  sediment  (when  the  dimension¬ 
less  grain  size  d*  =  [(s  —  llg/vW^o  >  15,  in  which  g  is  the 
acceleration  of  the  gravity  and  v  is  the  kinematic  viscosity  of 
water). 

Validation  of  the  Hypothesis 

To  validate  the  two  main  hypotheses  underlying  the  pro¬ 
posed  formula,  a  comparison  was  made  between  the  observed 
suspended  load  and  Equation  (4)  with  the  fitted  values  to  the 
observed  data  on  cR  and  e  for  each  experimental  case  (least 
squares  method  used  to  fit  the  exponential  profile).  For  the 
data  with  current  only  (see  Table  1),  very  good  agreement 
was  observed  for  all  data  sets.  In  Table  3,  the  main  results 
of  the  comparison  are  provided — i.e.,  the  percentage  of  data 
correctly  predicted  within  a  factor  of  2  or  5  (Pred.  X2  and 
Pred.  X5,  respectively)  and  the  mean  value  and  standard  de¬ 
viation  of  the  function  f(qj  =  log  |<3'sa,pred/«3ras,meaB I ,  where  gaa  pred 
and  QVmeas  are  the  predicted  and  measured  suspended  load 
values,  respectively.  The  table  shows  that  99%  of  the  data 
are  well  predicted  within  a  factor  of  2.  This  means  that  the 
assumptions  of  an  exponential  concentration  profile  and  a 
constant  velocity  over  depth  are  sufficiently  accurate  to  es¬ 
timate  the  suspended  load  in  the  steady  current  case. 

With  waves  and  current  combined,  the  results  are  not  as 
good  as  for  the  steady  current  data.  Only  66%  and  91%  of  the 
data  (72%  and  93%  in  the  case  of  breaking  waves)  are  pre¬ 
dicted  within  a  factor  of  2  and  5,  respectively  (see  Table  3). 
Indeed,  the  suspended  load  is  quite  often  underestimated. 
Large  uncertainties  occurred  when  estimating  the  suspended 
load,  both  for  the  actual  data  and  the  fitted  data.  Knowing 
the  concentration  profile  and  use  of  the  mean  current  velocity 
instead  of  the  velocity  profile  does  not  significantly  affect  the 
suspended  load.  It  induces  an  overestimation  of  less  than  5% 


for  a  logarithmic  profile  (longshore  current)  and  an  under¬ 
estimation  of  10%  to  an  overestimation  of  20%  in  more  com¬ 
plex  flows  (undertow).  The  most  important  assumption  for 
the  sediment  transport  seems  thus  to  be  the  exponential  con¬ 
centration  profile.  For  some  cases,  it  appears  that  a  power 
law  profile  fits  better  to  the  data,  especially  close  to  the  bed. 
For  many  cases — e.g.,  the  Grote  Speurwerk  (45m)  data  set — 
no  measurements  of  the  bed  features  are  available.  The  es¬ 
timation  of  the  measured  suspended  load  is  also  strongly 
linked  to  the  chosen  reference  level;  if  bed  forms  are  present, 
zR  is  considered  to  be  equal  to  half  of  the  measured  or  esti¬ 
mated  bed  form  height  (Van  Rijn,  1984b).  When  the  ripple 
characteristics  were  not  measured,  the  Soulsby  and  White- 
house  (2005)  formulas  were  used.  This  might  explain  the 
underestimation  observed  for  some  data  sets. 

SEDIMENT  DIFFUSIVITY 

The  sediment  diffusivity  e  is  a  fundamental  parameter  for 
estimating  the  concentration  profile.  It  is  a  function  of  bottom 
roughness  and  associated  shear  stress,  agitation  (mainly 
from  waves),  and  settling  velocity. 

A  General  Equation  for  the  Sediment  Diffusivity 

In  a  study  on  infilling  of  navigation  channels,  Kraus  and 
Larson  (2001)  employed  an  exponential  concentration  pro¬ 
file  to  estimate  the  suspended  load  transport  and  assumed, 
following  Battjes  (1983),  that  vertical  mixing  is  proportional 
to  the  energy  dissipation  of  wave  breaking.  Then,  to  employ 
a  general  formula  for  sediment  diffusivity,  it  seems  natural 
to  assume  that 


where  D  is  the  total  effective  dissipation  and 

D  =  k3bDb  +  k3cDc  +  klDw  (7) 

in  which  the  energy  dissipation  because  of  wave  breaking 
(Db)  and  bottom  friction  from  current  (Dc)  and  waves  (. Dw)  are 
simply  added  and  kb,  kc,  and  are  coefficients.  The  coeffi¬ 
cient  kb  mainly  corresponds  to  an  efficiency,  whereas  kc  and 
are  related  to  the  Schmidt  number  (ratio  between  the  ver¬ 
tical  eddy  diffusivity  of  the  particle  e  and  the  vertical  eddy 
viscosity  vv).  Typically,  Db  >  Dw  >  Dc,  and  in  many  cases, 
only  the  largest  dissipation  needs  to  be  considered.  Still,  the 
formula  for  e  is  simplified  because  mixing  should  vary 
through  the  water  column. 

The  energy  dissipation  in  the  bottom  boundary  layer  be¬ 
cause  of  a  current  can  be  written 

Dc  =  Tcw*c  (8) 

where  tc  and  w*c  are  bottom  shear  stress  and  shear  velocity 
from  the  current  only,  respectively.  This  expression  deviates 
from  the  standard  way  of  defining  dissipation  by  a  current, 
which  should  be  expressed  as  the  product  between  the  shear 
stress  tc  and  the  mean  velocity  Uc.  However,  the  use  of  m*c 
instead  of  Uc  yields  the  same  result  as  the  classical  mixing 
length  approach  (Equation  [9]). 
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Figure  3.  Estimation  of  the  Schmidt  number  ac  as  a  function  of  the  ratio 
WJu* c-  For  a  color  version  of  this  figure,  see  page  682. 


(D\m 

ec  =  kjr —  I  h  =  kcu%  ch  (9) 

Similarly,  the  instantaneous  energy  dissipation  in  the  bot¬ 
tom  boundary  layer  because  of  wave  motion  Z)wt(t)  can  be  ex¬ 
pressed  as  the  product  between  the  instantaneous  bed  shear 
stress  t wt(f)  and  the  instantaneous  shear  velocity  u*wt(f).  The 
sediment  diffusivity  as  the  result  of  waves  can  thus  be  writ¬ 
ten  as  an  average  over  the  wave  period: 

ew  =  —  j  h  =  k„u^„h  (10) 

where  Dw  (=  t wm*w),  tw,  and  u*w  are  the  wave-related  energy 
dissipation,  maximum  bottom  shear  stress,  and  maximum 
shear  velocity,  respectively. 

Effects  of  a  Steady  Current 

For  suspended  sediment  in  a  steady  current,  Rouse  (1938) 
observed  that  a  parabolic  equation  for  the  vertical  sediment 
diffusivity  best  fits  the  experimental  data.  Assuming  that  the 
parabolic  profile  is  a  correct  approximation  for  vertical  sedi¬ 
ment  diffusivity,  the  relationship  in  Equation  (11)  is  ob¬ 
tained, 


K  =  -JK  (11) 

6 

where  crc  is  the  Schmidt  number  (also  defined  as  a  |3-factor 
by  Nielsen  and  Teakle  [2004]  or  Van  Rijn  [1984b]).  In 
general,  uc  is  supposed  to  be  a  constant  and  ctc  =  1.  However, 
Sumer  and  Deigaard  (1981)  and  Van  Rijn  (1984b)  pointed 
out  that  centrifugal  forces  in  fluid  eddies  cause  sediment 
grains  to  be  thrown  outside  of  the  fluid  eddies,  which  increas¬ 
es  uc.  Another  reason  presented  by  Fredspe  and  Deigaard 
(1992,  pp.  231-234)  for  why  ctc  might  deviate  from  unity  is 


Table  4.  Prediction  (Pred.)  of  the  Schmidt  number  for  the  case  of  a  steady 
current  [^aj  =  log  \acpr,d/<Tcmea,\]. 


Authors 

Pred.  (%) 

f(cr  J 

X2 

X5 

Mean 

SD 

Van  Rijn  (1984b) 

78 

100 

-0.02 

0.25 

Rose  and  Thorne  (2001) 

75 

98 

-0.02 

0.32 

Equation  (12) 

76 

99 

-0.04 

0.26 

sediment  settling  out  of  the  surrounding  water  before  the  wa¬ 
ter  loses  its  earlier  composition  by  mixing.  Rose  and  Thorne 
(2001)  added  that  the  estimation  of  crc  might  be  affected  by 
the  settling  velocity,  which  varies  because  of  turbulence  (a 
smaller  particle  settles  faster  in  turbulence;  see  Murray, 
1970;  Nielsen,  1993).  Finally,  Nielsen  and  Teakle  (2004) 
proposed  a  Fickian  diffusivity  model  showing  the  effect  of  the 
size  of  the  particle  on  the  Schmidt  number.  They  also  argued 
that  the  Schmidt  number  might  be  less  than  unity  for  fine 
particles  over  a  flat  bed  because  the  mixing  length  could  be 
smaller  for  these  particles  than  for  the  fluid  (Muste  and  Pa¬ 
tel,  1997).  To  maintain  a  consistent  approach  between  the 
analysis  of  different  data  sets,  the  von  Karman  constant  was 
assumed  to  equal  the  clear  water  value  =  0.4,  even  though 
some  evidence  suggests  that  this  parameter  is  reduced  by  the 
presence  of  suspended  sediment  (on  the  basis  of  the  assump¬ 
tion  that  the  Schmidt  number  equals  1).  These  variations  are 
thus  included  in  the  Schmidt  number. 

On  the  basis  of  measurements  by  Coleman  (1981),  Van 
Rijn  ( 1984b)  suggested  an  expression  for  ctc  that  is  a  function 
of  the  ratio  between  the  settling  velocity  and  the  current- 
related  shear  velocity.  With  the  exponential  law  to  derive  the 
sediment  diffusivity,  it  is  also  possible  to  estimate  the  coef¬ 
ficient  ctc  from  the  experimental  concentration  profiles.  The 
experimental  results  are  presented  in  Figure  3  together  with 
the  Van  Rijn  and  Rose  and  Thorne  formulas  and  the  new 
equation  proposed  in  this  paper  (Equation  [12]),  as  well  as 
additional  data  from  Rose  and  Thorne  (2001)  and  Van  Rijn 
(1984b). 

To  put  forward  a  relationship  that  gives  physically  mean¬ 
ingful  results  for  all  cases,  it  should  be  considered  that  the 
Van  Rijn  equation  is  correct  only  for  WJu*c  <  1.  For  WJu *c 
»  1,  the  Schmidt  number  must  be  equal  to  unity  because 
suspension  is  negligible  (sediments  do  not  affect  the  flow). 
Thus,  a  new  expression  for  ac  is  proposed  in  Equation  (12) 
(see  also  Figure  3), 


'1“+ S) 

1  +  (Acl  +  Ac2  -  Dsin2^^j 


W*c 


1 


W*c 


>  1 


(12) 


where  Acl  =  0.40  andAc2  =  3.5. 

Overall,  the  new  equation  yields  good  predictions  compared 
with  the  data  (see  Table  4).  Nevertheless,  data  from  the  Peet 
(1999,  see  Van  Rijn  et  al.,  2001)  experiments  are,  in  general, 
overestimated,  perhaps  because  measurements  were  only 
carried  out  close  to  the  bed  (z  <  h/5). 
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Figure  4.  Estimation  of  the  Schmidt  number  crw  with  Equation  (14)  as 
a  function  of  the  ratio  Ws/u*w,  with  the  roughness  ratio  ka/d50  emphasized. 
For  a  color  version  of  this  figure,  see  page  682. 


Effects  of  Nonbreaking  Waves 

Following  the  results  obtained  for  the  steady  current,  the 
relationship  in  Equation  (13)  for  is  obtained  (the  factor  2/tt 
results  from  time  averaging  assuming  a  sinusoidal  wave), 


K  = 


(13) 


where  crw  is  the  wave-related  Schmidt  number. 

With  the  data  sets  presented  in  Table  2,  the  total  shear 
velocity  was  estimated  assuming  that  the  shear  velocities  due 
to  the  current  and  waves  can  be  added  linearly.  For  most  of 
the  cases  (wave  flumes),  the  shear  stress  from  the  current  is 
much  smaller  than  the  shear  stress  from  the  waves  and,  thus, 
can  be  neglected.  The  total  roughness  height  k„  was  estimat¬ 
ed  by  the  method  of  Soulsby  (1997),  adding  the  grain-relat¬ 
ed,  form  drag,  and  sediment  transport  components  ksg  =  2 d50, 
ks[,  and  kss,  respectively.  The  roughness  height  ks[  was  ob¬ 
tained  with  the  Kjm  (2004)  formula,  which  is  a  function  of 
the  bed  form  characteristics,  whereas  kBS  was  obtained  with 
the  Wilson  (1966,  1989)  formula. 

Following  the  method  used  for  sediment  diffusivity  because 
of  a  current  (see  previous  section),  the  correction  factor 
(Schmidt  number)  was  assumed  to  be  a  function  of  the  ratio 
WJu*„.  The  general  tendency  of  the  data  appears  similar  to 
the  steady  current  case  (see  Figures  3  and  4), 


CT„ 


(i T  W  \ 
Awi  +  Aw2sin2  —  - 

\2  «*„/ 

1  +  (Awl  +  Aw2  - 


if 


<  1 


if 


ws 


>  1 


(14) 


where  Awl  =  0.15  and  Aw2  =  1.5. 


Table  5.  Prediction  (Pred.)  of  the  diffusivity  for  waves  only  /f(e,J  = 
log 


Authors 

Pred.  (%) 

A* 

w) 

X2 

X5 

Mean 

SD 

Dally  and  Dean  (1984) 

14 

51 

0.67 

0.32 

Van  Rijn  (1993) 

38 

76 

0.29 

0.52 

Nielsen  (1992) 

56 

90 

-0.31 

0.57 

Equations  (10)  and  (14) 

69 

100 

0.09 

0.28 

The  wave-induced  Schmidt  number  (see  Equation  [14])  is 
often  much  smaller  than  that  found  for  a  steady  current. 
However,  because  the  friction  velocity  from  waves  is  gener¬ 
ally  much  larger  than  that  from  a  current,  the  mixing  attrib¬ 
utable  to  waves  is  also  much  larger.  In  Figure  4,  there  seems 
to  be  a  relationship  between  ctw  and  the  roughness  ratio 
kjd50.  Because  this  roughness  ratio  (and  the  total  shear 
stress)  is  calculated  from  empirical  formulas,  and  not  esti¬ 
mated  directly  from  the  experimental  data  (contrary  to  the 
data  with  current  only),  the  relationship  between  ew  and 
Ws/wH!w  exhibits  larger  scatter. 

The  results  obtained  by  Equations  (10)  and  (14)  are  better 
than  those  for  the  other  studied  formulas  (see  Table  5).  Thus, 
69%  of  the  data  is  correctly  predicted  within  a  factor  of  2  and 
100%  within  a  factor  of  5.  The  percentages  of  values  obtained 
within  a  factor  of  2  or  5,  as  well  as  the  mean  value  and  the 
standard  deviation  of  the  function  /Tew)  =  log|ewpred/ewmeas|, 
are  presented  in  Table  5. 

It  appears  that  of  the  existing  formula,  the  one  proposed 
by  Dally  and  Dean  (1984),  who  assumed  that  the  Rouse 
(1938)  expression  could  be  used,  presents  the  least  scatter, 
even  if  the  formula  in  general  overestimates  sediment  diffu¬ 
sivity.  This  expression  might  be  correct,  but  the  Schmidt 
number  appears  to  be  much  smaller  than  1.  The  more  com¬ 
plex  formulas  introduced  by  Nielsen  (1992,  pp.  215-217)  or 
Van  Rijn  (1993)  yield  better  results  (especially  for  the  Niel¬ 
sen  formula),  but  also  larger  scatter.  The  Nielsen  formula 
tends,  however,  to  largely  underestimate  the  results  when 
UJWS  >  18. 

Wave-Current  Interaction 

Simply  adding  the  sediment  diffusivity  from  the  current 
and  waves  (Equations  [9],  [10],  [12],  and  [14])  leads  to  over¬ 
estimation  compared  with  the  data.  This  overestimation 
could  be  because  the  Schmidt  number  should  be  the  same  for 
the  current  as  for  the  waves.  A  more  physical  description  of 
the  wave  and  current  interaction  should  be  based  on  a  unique 
Schmidt  number,  calculated  as  an  empirical  weighted  value 
between  ctc  and  aw, 

Hew  =  +  (1  ~  ^t)frw  (15) 

where  Xt  =  0C/(0C  +  0W),  in  which  0C  =  t,/[(s  —  1  )gd50],  and  0W 
=  tw/[(s  -  l)gii50],  which  are  the  current-related  and  wave- 
related  Shields  parameters,  respectively. 

Table  6  shows  the  results  obtained  with  the  studied  for¬ 
mulas.  Because  larger  values  are  observed  in  cases  with  wave 
and  current  interaction,  the  formulas  in  which  the  sediment 
diffusivity  was  overestimated  for  waves  only  provide  better 
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Table  6.  Prediction  (Pred.)  of  the  sediment  diffusivity  in  the  wave  and 
current  interaction  case  [i (ecJ  =  log  I  eCUIpTedl  ecwmeas  \J. 


Pred.  (%) 

f(£c 

,) 

Authors 

X2 

X5 

Mean 

SD 

Dally  and  Dean  (1984) 

47 

88 

0.26 

0.39 

Van  Rijn  (1993) 

38 

83 

0.19 

0.50 

Nielsen  (1992) 

00 

01 

-1.68 

0.80 

Equation  (9)  +  Equation  (10) 
idem  with  ac  =  crw  =  acw 

50 

85 

0.31 

0.36 

[Equation  (15)1 

65 

92 

0.08 

0.39 

results  (see  formulas  by  Dally  and  Dean,  1984;  Van  Rijn, 
1993).  The  Nielsen  (1992)  formula  yields  poor  results  be¬ 
cause  UJWS  >  18  for  most  of  the  cases.  The  proposed  formula 
to  calculate  ecw  by  Equation  (15)  yields  the  best  results,  al¬ 
though  it  often  overestimates  when  large  bed  roughness  is 
computed. 

Effects  of  Breaking  Waves 

In  the  case  of  breaking  waves,  the  energy  dissipation  was 
calculated  with  the  energy  dissipation  of  a  bore  analogy 
(Svendsen,  1984)  with  the  coefficient  Ae  =  2  tanh(5£«,)  pro¬ 
posed  by  Stive  (1984)  to  take  into  account  breaker-type  ef¬ 
fects  (t  =  m/VHjIZ  is  the  Irribaren  parameter  defined  for 
deep  water,  m  is  the  mean  slope  of  the  beach,  and  H„  and  Lx 
are  the  deep-water  wave  height  and  length,  respectively).  In 
the  random  waves  case,  the  coefficient  ab  =  exp[— (-ybA/H'rmB)2] 
should  be  added  to  take  into  account  the  percentage  of  break¬ 
ing  waves  (see  Larson,  1995),  where  yb  is  the  breaker  depth 
index  and  Hrms  is  the  root  mean  square  wave  height,  neglect¬ 
ing  breaking.  As  a  first  approximation,  the  efficiency  coeffi¬ 
cient  was  found  to  be  constant:  kb  =  0.010.  Even  if  some  dis¬ 
persion  exists  for  the  compiled  data,  Equation  (6)  yields  re¬ 
sults  as  good  as  for  the  nonbreaking  cases:  72%  (96%)  of  the 
data  are  well  predicted  within  a  factor  of  2  (5),  and  the  mean 
and  standard  deviation  of  the  function  f(eb)  =  log|eb  pred/eb  meaa| 
is  0.03  and  0.32,  respectively. 

REFERENCE  CONCENTRATION 

Einstein  (1950)  proposed  that  the  reference  concentration 
could  be  a  function  of  the  bed  load  transport.  Following  Mad¬ 
sen,  Tajima,  and  Ebersole  (2003),  the  reference  volumetric 
bed  concentration  can  be  estimated  from  the  volumetric  bed 
load,  assuming  qs  =  cR  U„  where  U„  is  the  velocity  of  the  bed 
load  layer.  The  bed  load  can  be  written  following  the  results 
of  Camenen  and  Larson  (2005),  namely  qs  <*  03/2  exp(-4.50/ 
0cr),  where  0  is  the  Shields  parameter  (the  subscript  “cr”  in¬ 
dicates  its  critical  value  for  the  inception  of  sediment  motion; 
see  Camenen  and  Larson,  2005;  Soulsby,  1997).  Madsen, 
Tajima,  and  Ebersole  (2003)  proposed,  as  a  first  approxi¬ 
mation,  that  the  speed  of  the  bed  load  layer  is  proportional 
to  the  shear  velocity,  UB  «  0V2.  The  bed  reference  concentra¬ 
tion  could  thus  be  written  as  Equation  (16), 

-4.5^  (16) 

“cr  / 


Table  7.  Prediction  (Pred.)  of  the  reference  concentration  assuming  an 
exponential  sediment  concentration  profile  for  the  studied  data  set  with 
steady  current  only  [f(cB)  =  log(cRp„J cRmepJ], 


Pred.  (%) 

f(c  r) 

Authors 

X2 

X5 

Mean 

SD 

Madsen,  Tajima,  and 
Ebersole  (2003) 

27 

50 

0.75 

0.83 

Nielsen  (1986,  1992) 

13 

50 

0.64 

0.43 

Equations  (16)  and  (17) 

38 

83 

-0.12 

0.55 

Equation  (16)  with 

Acr  =  5  X  10~4 

28 

65 

0.05 

0.77 

where  0T  is  the  transport-dependent  Shields  parameter  and 
0M  is  the  maximum  Shields  parameter.  For  current  only,  0M 
=  0T  =  0C. 

Effect  of  a  Current 

For  the  data  with  a  steady  current,  the  coefficient  AcE  was 
found  to  vary  from  5  X  10~6  to  4  X  10~2.  The  use  of  a  constant 
mean  value  AcR  =  5  X  1CL4  produces  results  already  better 
than  the  Madsen,  Tajima,  and  Ebersole  (2003)  formula 
(see  Table  7).  Van  Rijn  (1984a,  1984b)  observed  that  the  ref¬ 
erence  concentration  cE  is  a  function  of  the  dimensionless 
grain  size  d*,  but  to  a  varying  power  depending  on  the  pres¬ 
ence  or  absence  of  bed  forms.  For  the  compiled  data  (see  Ta¬ 
ble  1),  the  dimensionless  grain  size  d*  varies  from  1  to  18. 
Improved  results  were  obtained  by  calibrating  AcR  as  a  func¬ 
tion  of  the  dimensionless  grain  size,  as  in  Equation  (17). 

AcR  =  1.5  X  10-3  exp(-0.2d*)  (17) 

The  percentages  for  values  obtained  within  a  factor  of  2  or 
5,  as  well  as  the  mean  value  and  the  standard  deviation  of 
the  ratio  f(cR)  =  log(cRpred/cEmeaa),  are  presented  in  Table  7. 
It  appears  that  the  Madsen,  Tajima,  and  Ebersole  (2003) 
formula  presents  correct  results  compared  with  the  experi¬ 
mental  values  of  cE  assuming  an  exponential  profile.  How¬ 
ever,  this  formula  seems  not  to  be  sensitive  enough  to  the 
Shields  parameter:  it  gives  more  or  less  a  constant  value  for 
each  data  set.  The  Nielsen  (1986,  1992,  pp.  201-233)  for¬ 
mula,  although  fitted  with  data  on  waves  only,  presents  cor¬ 
rect  results  for  the  laboratory  experiments  (data  from  Bar¬ 
ton  and  Lin,  1955;  Laursen,  1958;  Peet,  1999  [see  Van 
Rijn  et  al.,  2001]).  It  tends,  however,  to  overestimate  the  re¬ 
sults  for  the  field  experiments,  even  if  bed  form  influence  on 
the  Shields  parameter  were  not  taken  into  account  because 
information  was  not  available. 

In  Figure  5,  the  predicted  reference  concentration  cE  from 
Equations  (16)  and  (17)  is  plotted  against  the  estimated  ref¬ 
erence  concentration  assuming  an  exponential  profile.  Even 
if  the  results  are  in  agreement  overall,  it  seems  that  sensi¬ 
tivity  to  the  Shields  parameter  should  be  larger.  The  refer¬ 
ence  concentration  is  generally  overestimated  for  low  shear 
stresses.  The  prediction  of  the  reference  concentration  is  sig¬ 
nificantly  improved  compared  with  the  Madsen,  Tajima, 
and  Ebersole  (2003)  formula:  nearly  40%  of  the  data  are 
predicted  within  a  factor  of  2  and  85%  within  a  factor  of  5. 
The  mean  value  of  f(cR)  is  closer  to  zero,  and  its  standard 


Cr  =  AcR0Texp 
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Figure  5.  Predicted  reference  concentration  cR  with  Equations  (16)  and 
(17)  vs.  the  estimated  reference  concentration  assuming  an  exponential 
profile.  For  a  color  version  of  this  figure,  see  page  682. 


deviation  is  reduced  compared  with  the  previous  formulas. 
However,  some  dispersion  appears  for  the  data  sets  of  Peet 
(1999,  see  Van  Rijn  et  al.,  2001)  and  Voogt  et  al.  (1991), 
which  produces  the  negative  value  for  the  mean  of  f(cR). 

Effect  of  Nonbreaking  Waves 

For  the  waves  only  case,  according  to  the  results  of  Ca- 
menen  and  Larson  (2005)  for  bed  load  transport,  the  mean 
shear  stress  0w  m  is  used  for  the  transport-dependent  term 
0T  =  9w,m>  whereas  the  maximum  wave  shear  stress  0W  is 
used  for  the  critical  shear  stress  (0M  =  0W).  The  mean 
Shields  parameter  attributable  to  the  waves  is  defined  as 
0w,m  =  %/"„/[(s  -  l)g'<f50][fow  ww(f)2  dt],  whereas  the  maximum 
Shields  parameter  is  defined  as  0W  =  Vif^UiUis  —  l)gG?50], 
where  /w  is  the  wave-related  friction  coefficient  and  0m  w  = 
O.50w  in  the  case  of  a  sinusoidal  wave. 

For  the  data  with  waves  prevailing  (|t/c|  <  0.05  m/s;  see 
Table  2),  a  comparison  was  made  between  the  different  for¬ 
mulas  studied  and  the  data  (see  Table  8  and  Figure  6).  Equa¬ 
tion  (16)  with  Equation  (17)  for  current  only  still  presents 
correct  results,  although  the  dispersion  is  larger.  The  effect 
of  grain  size  seems  not  to  be  as  significant  as  for  the  results 
with  current  only.  Similar  results  are  obtained  with  Equation 
(16)  in  which  AcE  =  5  X  1CL4.  However,  the  range  of  values 
on  d *  was  larger  for  the  data  set  with  current  only  ( d *  from 
1.0  to  18).  In  the  case  of  current  only,  d*  <  5  for  40%  of  the 
data,  whereas  95%  of  the  data  was  below  this  value  for  waves 
only.  This  result  could  explain  the  difference  in  the  results 
for  the  current,  and  especially  the  difference  observed  with 
the  use  of  Equation  (16)  and  AcE  =  5  X  1CL4  (overestimation 
for  the  current  data  set  and  underestimation  for  the  wave 
data  set). 

Another  reason  for  the  large  scatter  comes  from  the  uncer¬ 
tainties  in  the  bed  form  characteristics  (measured  or  esti¬ 
mated)  and  in  the  total  shear  stress  calculation.  AcE  is  ob- 


Table  8.  Prediction  (Pred.)  of  the  reference  concentration  with  the  waves- 
only  data  set  [\\]c\  <  0.05  m/s,  i(CR)  =  log(cRpredlcRml!as)]. 


Pres.  (%) 

f(c  r) 

Authors 

X2 

X5 

Mean 

SD 

Nielsen  (1986,  1992) 
Madsen,  Tajima,  and 

26 

60 

0.34 

1.04 

Ebersole  (2003) 

20 

45 

0.76 

0.57 

Equations  (16)  and  (17) 
Equation  (16)  with 

31 

66 

-0.27 

0.63 

AIr  =  5  X  1Q-4 

31 

67 

-0.26 

0.61 

served  to  be  a  function  of  the  ripple  height,  Ht,  or,  more  spe¬ 
cifically,  of  the  roughness  height  ratio  kjdm  (see  Figure  6). 
Finally,  the  Madsen,  Tajima,  and  Ebersole  (2003)  formula 
(as  well  as  Equation  [16]  with  AcE  =  5  X  10  4)  again  shows 
reasonable  results  because  it  is  not  as  sensitive  to  d *  and  0. 
In  Figure  6,  the  roughness  height  ratio  is  emphasized  and 
shows  that  the  calculation  of  the  total  shear  stress  induces 
large  uncertainties  (assuming  that  the  reference  concentra¬ 
tion  cE  should  not  be  a  function  of  the  ripple  height  or  the 
roughness  ratio,  but  only  of  the  total  shear  stress).  It  appears 
that  the  more  overestimated/underestimated  the  reference 
concentration  is,  the  higher/smaller  the  roughness  ratio. 

Wave-Current  Interaction 

For  the  wave-current  interaction,  the  intuitive  Shields  pa¬ 
rameters  to  be  used  in  Equation  (16)  are  0T  =  0cwm  and  0M 
=  0CW.  To  simplify  the  calculations,  the  mean  and  maximum 
Shields  parameters  from  the  wave-current  interaction  can  be 
obtained  simply  through  addition:  i.e.,  0cw  m  =  (02  +  02  m  + 
20w  m0c  cos  tp)4'2  and  0CW  =  (02  +  02  +  20W0C  cos  9 )y2,  respec¬ 
tively,  where  9  is  the  angle  between  the  wave  and  current 
directions. 


Figure  6.  Reference  concentration  cR  estimated  from  the  data  compiled 
(with  waves  only,  \UC\  <  0.05  m/s)  vs.  cR  calculated  with  Equations  (16) 
and  (17)  (roughness  ratio  emphasized).  For  a  color  version  of  this  figure, 
see  page  682. 
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Table  9.  Prediction  (Pred.)  of  the  reference  concentration  from  the  studied 
data  set  with  wave-current  interaction  [parentheses  hold  the  results  for 
breaking  waves  only.  f(cR)  =  log(cRpredl cRmraJJ. 


Pred.  (%) 

f(c r) 

Authors 

X2 

X5 

Mean 

SD 

Nielsen  (1986,  1992) 
Madsen,  Tajima,  and 

42  (30) 

68  (60) 

-0.25  (0.30) 

0.73  (0.90) 

Ebersole  (2003) 

02  (08) 

10  (37) 

1.11  (0.90) 

0.43  (0.50) 

Equations  (16)  and  (17)  52  (62) 
Equations  (16)  with 

82  (90) 

0.17  (0.01) 

0.53  (0.43) 

Alr  =  5  x  10-4 

57  (59) 

85  (88) 

0.07  (-0.04) 

0.50  (0.45) 

The  percentages  for  values  obtained  within  a  factor  of  2  or 
5,  as  well  as  the  mean  value  and  the  standard  deviation  of 
the  ratio  f(cR)  =  log(cRpred/cRexp),  for  the  wave-current  inter¬ 
action  are  presented  in  Table  9  for  the  investigated  formulas 
and  Equation  (16).  The  Nielsen  formula  presents  better  re¬ 
sults  compared  with  the  waves  only  case,  but  it  is  still  very 
scattered  (because  it  is  a  function  of  the  ripple  characteris¬ 
tics).  On  the  other  hand,  the  Madsen,  Tajima,  and  Eber- 
sole  (2003)  formula  does  not  have  so  much  scatter,  but  it 
generally  overestimates,  as  does  the  present  formula  with  0T 
=  0CW  m.  However,  the  new  formula  still  presents  the  best  re¬ 
sults  among  those  studied.  It  appears  that  the  computation 
of  the  mean  wave  and  current  Shields  parameter  0cw  m  sig¬ 
nificantly  influences  the  results.  Compared  with  the  waves 
only  cases,  the  formula  yields  better  results  but  generally 
overestimates.  The  use  of  Equation  (16)  with  AcR  =  5  X  ICh4 
produces  surprisingly  similar  estimates. 

Effect  of  Breaking  Waves 

As  a  first  approach,  it  can  be  assumed  that  wave  breaking 
does  not  affect  the  reference  concentration,  but  only  the  sed¬ 
iment  diffusivity.  Indeed,  as  shown  by  Nielsen  (1992,  p. 
219),  the  turbulence  induced  by  the  breakers  generally  occurs 
in  the  upper  part  of  the  water  column;  it  should  not  influence 
the  bottom  concentration  significantly.  The  reference  concen¬ 
tration  could,  however,  be  enhanced  by  the  breakers  in  the 
case  of  plunging  waves,  in  which  the  generated  turbulent  jet 
might  penetrate  to  the  bottom. 

The  data  sets  presented  in  Table  2  involve  many  experi¬ 
mental  cases  in  which  breaking  waves  occurred.  Table  9  pre¬ 
sents  the  prediction  results  depending  on  the  chosen  formula. 
Even  if  the  results  are  scattered,  Equations  (16)  and  (17)  pre¬ 
sent  the  best  results  among  the  formulas  studied,  with  62% 
of  the  data  correctly  predicted  within  a  factor  of  2  and  90% 
within  a  factor  of  5.  No  clear  effect  from  the  type  of  breaker 
could  be  observed  from  the  collected  data. 

SUSPENDED  LOAD 

The  use  of  Equation  (4)  together  with  the  expressions  for 
sediment  diffusivity  (Equation  [6])  and  the  reference  concen¬ 
tration  (Equation  [16])  allows  for  prediction  of  the  suspended 
load  for  a  steady  current,  waves  and  current  combined,  and 
breaking  waves. 


Table  10.  Prediction  (Pred.)  of  suspended  load  transport  in  the  steady 
current  case  /ifq„J  =  log(q_SBp„dl  q_ssmeas)]. 


Authors 

Pred.  (%) 

f(qj 

X2 

X5 

Mean 

SD 

Bijker  (1968) 

24 

45 

0.60 

1.04 

Engelund  and  Hansen  (1972) 

31 

55 

0.65 

0.85 

Bailard  (1981) 

33 

72 

0.32 

0.69 

Van  Rijn  (1984b) 

30 

69 

-0.27 

0.98 

This  work  with  equation  (16) 

37 

79 

-0.10 

0.57 

Comparison  with  Data  for  Current  Only 

A  comparison  between  the  predicted  and  observed  sus¬ 
pended  sediment  load  is  presented  in  Table  10  and  Figure  7. 
In  general,  the  proposed  formula  (Equation  [4])  shows  correct 
behavior.  The  obtained  results,  however,  appear  to  be  highly 
dependent  on  the  estimate  of  the  reference  concentration 
when  a  large  dispersion  is  observed.  Indeed,  the  equation  pro¬ 
posed  for  cR  for  a  steady  current  only  (i.e.,  Equations  [1]  and 
[17])  produces  much  better  results  than  the  use  of  a  constant 
value  for  AcR  (=  5  X  10  4).  An  increase  in  accuracy  of  nearly 
10%  and  a  decrease  in  the  standard  deviation  by  10%  can  be 
observed. 

In  Figure  7,  similar  behavior  as  that  of  the  reference  con¬ 
centration  predictions  is  observed:  a  general  overestimation 
for  the  Anderson  (1942),  Peet  (1999,  see  Van  Rijn  et  al., 
2001),  and  Scott  and  Stephens  (1966)  data  sets  and  a  gen¬ 
eral  underestimation  for  the  Barton  and  Lin  (1955)  and 
Laursen  (1958)  data.  A  comparison  with  other  semiempiri- 
cal  formulas  found  in  the  literature  showed  that  the  proposed 
relationship  significantly  improves  the  results  when  Equa¬ 
tion  (16)  is  used  for  the  reference  concentration. 


Figure  7.  Comparison  between  the  observed  suspended  sediment  load 
and  the  calculated  suspended  sediment  load  with  Equation  (4)  and  the 
predicted  values  for  cE  (Equation  [16])  and  e  (Equation  [7])  for  the  current 
only  case.  For  a  color  version  of  this  figure,  see  page  682. 
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Table  11.  Prediction  (Pred.)  of  suspended  load  transport  in  the  interac¬ 
tion  between  current  and  waves  [parentheses  show  the  results  for  breaking 
waves  only,  f(q „)  =  log( <i,s.pred 

!  Qss.meosH 


Pred.  (%) 

f(qj 

Authors 

X2 

X5 

Mean 

SD 

Bijker  (1968) 

10  (23) 

40  (59) 

0.83  (0.43) 

0.60  (0.70) 

Bailard  (1981) 

19  (30) 

58  (74) 

0.65  (0.47) 

0.52  (0.51) 

Van  Rijn  (1993) 

30  (23) 

70  (62) 

-0.22  (-0.03) 

0.74  (0.69) 

Equations  (4),  (6), 
and  (16) 

44  (42) 

77  (76) 

0.23  (-0.05) 

0.61  (0.60) 

Comparison  with  Data  for  Waves  and  Current 
Combined 

A  comparison  between  the  predicted  and  observed  sus¬ 
pended  sediment  load  for  the  wave-current  interaction  is  pre¬ 
sented  in  Table  11  and  Figure  8.  It  appears  that  the  proposed 
formula  (Equation  [4])  presents  overall  good  results.  Again, 
the  results  are  highly  dependent  on  the  estimation  of  the  ref¬ 


erence  concentration  and,  thus,  as  shown  in  the  Reference 
Concentration  section,  on  the  estimation  of  the  roughness 
height  and  total  shear  stress.  Moreover,  plunging  breaking 
waves  can  induce  larger  values  on  the  reference  concentra¬ 
tion,  as  discussed  in  the  Effect  of  Breaking  Waves  section.  Be¬ 
cause  the  prediction  of  sediment  diffusivity  is  generally  less 
scattered,  if  an  overestimation/underestimation  is  observed 
for  the  prediction  of  the  reference  concentration  (Grote 
Speurwerk  (35m),  Vessem,  Bayram  et  al.  [2001]/Kroon 
[1991]  data  sets),  the  same  observation  can  be  made  for  the 
resulting  suspended  load. 

The  percentages  for  values  obtained  within  a  factor  of  2 
and  5,  as  well  as  the  mean  value  and  the  standard  deviation 
of  the  ratio  f(qj  =  log(qs  pred/qs  exp),  are  presented  in  Table  11 
for  nonbreaking  and  breaking  waves  (results  in  parentheses 
are  for  breaking  waves).  For  the  wave-current  interaction 
without  breaking  waves,  the  Van  Rijn  (1993)  formula  yield 
the  best,  though  underestimated,  results.  The  large  number 
of  parameters  and  their  relative  complexity  might  explain  the 
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Figure  8.  Comparison  between  observed  and  predicted  values  of  the  suspended  sediment  load  in  the  wave-current  interaction  case  according  to  Equation 
(4).  For  a  color  version  of  this  figure,  see  page  682. 
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observed  scatter  because  it  is  more  sensitive  to  any  given 
parameter.  The  result  obtained  with  the  Van  Rijn  (1993) 
formula  is  poorer  when  waves  are  breaking.  On  the  other 
hand,  it  appears  that  the  Bailard  (1981)  formula  presents 
the  best  results  among  the  studied  formulas  for  breaking 
waves.  This  formula  was  calibrated  to  estimate  the  suspend¬ 
ed  load  in  the  surf  zone,  thus  yielding  good  predictions  under 
breaking  waves.  The  Bailard  (1982)  formula  also  produces 
less  dispersion.  Indeed,  this  formula  is  not  sensitive  to  shear 
stress  (only  an  average  friction  coefficient  is  introduced),  and 
it  is  simple  enough  to  reduce  the  dispersion  of  the  results. 
However,  no  improvement  of  the  results  could  be  obtained 
because  the  formula  is  basically  only  a  function  of  the  current 
and  wave  velocities  at  the  bottom.  Calculation  of  the  friction 
coefficient  also  appeared  to  affect  the  results  of  the  Bailard 
(1982)  formula.  In  the  same  way  as  the  Bijker  (1968)  for¬ 
mula,  the  Bailard  (1982)  formula  generally  tends  to  over¬ 
estimate  the  prediction. 

Equation  (4)  yields  the  best  results  for  both  nonbreaking 
and  breaking  cases.  The  result  appears,  however,  to  overes¬ 
timate  for  the  case  of  nonbreaking  waves  and  to  underesti¬ 
mate  slightly  for  the  case  of  breaking  waves  (see  the  Kroon 
[1991]  data  set  in  Figure  8).  Improvements  in  the  results 
could  be  obtained  with  better  predictions  of  the  total  shear 
stress  and  the  reference  concentration,  as  shown  in  the  Ef¬ 
fects  of  Nonbreaking  Waves  section. 

CONCLUSION 

In  this  paper,  we  presented  a  semi-empirical  formula  to 
estimate  the  suspended  load  resulting  from  a  mean  current 
and  waves.  Comparison  between  formula  predictions  and 
measured  suspended  load  from  a  large  data  set  demonstrated 
that  the  formula  overall  produces  satisfactory  predictions  of 
the  transport  rate  for  a  wide  range  of  hydrodynamic  and  sed¬ 
iment  conditions.  Assuming  an  exponential  profile  for  the 
sediment  concentration  and  a  constant  value  over  depth  for 
the  time-averaged  current  velocity,  the  resulting  sediment 
transport  rate  can  be  estimated  from  a  simple  equation  (see 
Equation  [4]).  The  two  main  parameters  are  mean  sediment 
diffusivity  over  depth  and  the  bottom  reference  concentra¬ 
tion. 

A  relationship  for  sediment  diffusivity  is  proposed  assum¬ 
ing  a  linear  combination  of  mixing  by  breaking  waves  and 
energy  dissipation  in  the  bottom  boundary  layer  from  a  mean 
current,  waves,  or  both  (see  Equations  [6]  and  [7]).  In  the 
boundary  layer,  dissipation  from  the  current/waves  was  ex¬ 
pressed  as  the  product  between  a  force  (bottom  shear  stress) 
and  a  velocity  (shear  velocity)  so  that  it  is  consistent  with  the 
classical  mixing  length  approach.  Formulas  to  predict  the 
Schmidt  number  a  were  developed  for  a  mean  current  and 
waves  separately.  For  mixing  by  breaking  waves,  an  efficien¬ 
cy  coefficient  was  introduced,  and  its  value  was  determined 
through  calibration  with  experimental  data. 

Following  the  Madsen,  Tajima,  and  Ebersole  (2003)  ap¬ 
proach,  the  reference  concentration  was  found  to  be  propor¬ 
tional  to  the  mean  Shields  parameter,  including  the  effect  of 
the  critical  Shields  parameter  introduced  by  Camenen  and 
Larson  (2005).  The  results  displayed  considerable  scatter, 


mainly  because  of  uncertainties  in  the  prediction  of  the  total 
Shields  parameter,  especially  when  bed  forms  were  present. 
Also,  as  Van  Rijn  (1993)  suggested,  the  dimensionless  grain 
size  d*  was  taken  into  account  in  the  calculation  of  the  ref¬ 
erence  concentration.  It  significantly  improves  the  results  for 
cases  with  current  only  or  with  waves  only. 

The  resulting  formula  for  the  suspended  load  is  robust  and 
effective,  and  it  gives  the  best  results  among  the  formulas 
studied  in  most  cases,  although  some  dispersion  still  exists. 
Furthermore,  because  it  is  a  physically  based  formula,  an  im¬ 
provement  of  the  process  knowledge  (e.g.,  for  the  estimation 
of  the  total  shear  stress)  could  easily  be  taken  into  account 
in  the  formula  and  thus  improve  the  results.  When  the  wave- 
related  suspended  load  prevails  (phase  lag  effects  over  a  rip¬ 
pled  bed),  Equation  (5)  can  be  modified  to  take  into  account 
these  effects  by  adding  a  coefficient  upl  s  that  will  decrease 
the  characteristic  onshore  velocity  and  increase  the  charac¬ 
teristic  offshore  velocity. 
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